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A generalized description of Talbot-Lau interference with matter waves is presented, which ac- 
counts for arbitrary grating interactions and realistic beam characteristics. The dispersion interac- 
tion between the beam particles and the optical elements strongly influences the interference pattern 
in this near-field effect, and it is known to dominate the fringe visibility if increasingly massive and 
complex particles are used. We provide a general description of the grating interaction process by 
combining semiclassical scattering theory with a phase space formulation. It serves to systematically 
improve the eikonal approximation used so far, and to assess its regime of validity. 
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I. INTRODUCTION 



The ability of material particles to show wave-like in- 
terference is one of the central predictions of quantum 
mechanics. While the early experimental tests worked 
with elementary particles [l|, H] , interferometry of atoms 
is by now a matured field of physics [3, |3j la] ■ It is in par- 
ticular the ability to cool and to control atoms using laser 
techniques that has propelled atom interferometry into a 
versatile tool for precision measurements, e.g. [g,[3,[aQ. 

As far as more complex and more massive objects are 
concerned, we are currently witnessing this transition 
from the proof-of-principle demonstration of their wave 
nature to the use of interferometry for quantitative mea- 
surements. Specifically, the static [lO| and the dynamic 
[ll| bulk polarizability of fuUerene molecules was mea- 
sured recently with a molecule interferometer. Also the 
controlled observation of decoherence, due to collisions 
with gas particles [Ij] or due to the emission of thermal 
radiation [ij|, has been used to characterize the interac- 
tion strength (or cross-section) of fuUerenes with exter- 
nal degrees of freedom. Another motivation for study- 
ing interference with large molecules is to test quantum 
mechanics in unprecedented regimes by establishing the 
wave nature of ever more massive objects |14| . 

The above-mentioned interference experiments with 
large molecules are based on the near-field Talbot-Lau 
effect, where three gratings are used which serve, in turn, 
to produce coherence in the beam, to bring it to in- 
terference, and to resolve the fringe pattern. This is 
an established technique in atom and electron interfer- 
ometry [H, [H, m [il, [H, and it is the method of 
choice for massive and bulky molecules (such as fuUerenes 
[20|, meso-tetraphenylprorphyrins [2l|, or functionalized 
azobenzenes [22] )■ The main reason is that a Talbot- 
Lau interferometer (TLI) tolerates beams which are rel- 
atively weakly coUimated, thus alleviating the increasing 
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difficulty in producing brilliant beams if the particles get 
more complex. Another important advantage compared 
to far-field setups, which is essential if one wants to in- 
crease the particle mass by several orders of magnitude, 
is the favorable scaling behavior of Talbot-Lau interfer- 
ometers with respect to the de Broglie wavelength Q. 

As a specific feature of Talbot-Lau interferometers, the 
forces between the particles in the beam and the diffrac- 
tion grating influence the interference pattern much 
stronger than in a far fleld setup. This is due to the fact 
that the different diffraction orders do not get spatially 
separated in a TLL Rather, all the orders interfere among 
each other, producing a resonant recurrence of the pat- 
tern whenever the so-called Talbot condition is met. As 
a consequence, even tiny distortions of the matter wave 
may lead to significant changes of the fringe visibility — 
requiring, for example, the modification of the van der 
Waals force due to retardation effects [2^ to be taken 
into account when describing the diffraction of fuUerenes 
at gold gratings. 

This effect of the dispersion force between the parti- 
cle and the grating wall gets more important as the mass 
and structure of the molecule grows. In particular, it sets 
increasingly strict requirements on the monochromatic- 
ity, i.e, the velocity spread permissible in the molecular 
beam. A very recent development, undertaken to reduce 
this influence, is the Kapitza-Dirac Talbot-Lau Interfer- 
ometer (KDTLI), where the second material grating is 
replaced by thepure phase grating produced by a stand- 
ing light field |2a |. 

The available theoretical descriptions of molecular 
Talbot-Lau interference account for the grating forces in 
terms of a simple eikonal phase shift [2J, [2a, [2Q| (an 
expression originally derived by nuclear physicists as an 
asymptotic high-energy approximation to the Lippmann- 
Schwinger equation in scattering theory [23, |28]). This 
approximation ceases to be valid with a growing influ- 
ence of the particle-grating interaction, and, due to its 
non-perturbative nature, its range of validity is not easy 
to assess. Therefore, given the quest for testing quan- 
tum mechanics with ever larger particles and given the 
increased precision required in metrological applications. 



there is a clear need to extend the theoretical description 
of near-field matter wave interference beyond the eikonal 
approximation. 

The main purpose of this paper is therefore to develop 
a generalized formulation of the coherent Talbot-Lau ef- 
fect. By combining a scattering theory formulation with 
semiclassical approximations we incorporate the effect of 
the grating interaction systematically beyond the eikonal 
approximation. As a prerequisite for its implementation, 
the established theory of near field interference first needs 
to be extended to account for the effects of finite angu- 
lar dispersion in the molecular beam. We show how this 
can be done transparently by using the phase space for- 
mulation of quantum mechanics. As a by-product, this 
formulation permits us to quantify the adjustment preci- 
sion required in realistic experiments. 

The structure of the article is as follows. In Section HIl 
we develop a generalized theory of Talbot-Lau interfer- 
ence, which is formulated independently of the particular 
choice of how to incorporate the grating interaction. We 
also establish the relation to the previous treatments by 
applying the eikonal approximation. In Section IIIIl we 
review relevant realistic descriptions of the grating in- 
teraction, and numerically illustrate their effect in the 
eikonal approximation. Section HVl is devoted to the de- 
velopment of a semiclassical formalism to go beyond the 
elementary eikonal approximation. Its predictions are 
numerically evaluated and compared in Section |Vl using 
the experimental parameters of the molecular interfer- 
ence experiments carried out in Vienna [20|, |22| . Finally, 
we present our conclusions in Section [VII 



II. PHASE SPACE DESCRIPTION OF THE 
TALBOT-LAU EFFECT 

The effect of Talbot-Lau interference can be described 
in a particularly transparent and accessible fashion by 
using the phase space representation of quantum me- 
chanics [2^, [20, [31I, ^M- This is demonstrated in [2a |. 
where both the coherent effect and the consequences of 
environmental interactions are formulated in terms of the 
Wigner function of the matter wave beam. As also shown 
there, the analogy of the Wigner function with the clas- 
sical phase space distribution allows one to evaluate the 
predictions of classical and quantum mechanics in the 
same framework, a necessary step if one wants to distin- 
guish unambiguously quantum interference from a possi- 
ble classical shadow effect. 

However, the treatment in [26[ is based on a number 
of idealizations, which must be reconsidered in view of 
a more refined description of the particle-grating inter- 
action. The least problematic approximation is to dis- 
regard the motion in the direction parallel to the grat- 
ing slits, which is permissible due to the translational 
symmetry of the setup in this direction. It follows that 
a two-dimensional description involving the longitudinal 
motion (denoted by z) and the transverse motion (de- 



noted by x) of the beam is required in principle. In front 
of the interferometer these two degrees of freedom are 
well approximated by a separable state involving trans- 
verse momenta \p\ <^ Pz- If the grating interaction is 
treated in eikonal approximation, as done in [26|, this 
implies that the transverse and the longitudinal motion 
remain separable throughout. One can then resort to 
an effectively one-dimensional description, characterized 
by a fixed longitudinal momentum pz — h/X. The z 
coordinate then represents a time t = mz/pz and the 
longitudinal propagation of the beam along z effectively 
evolves the one-dimensional transverse beam state during 
the time t. The finite distribution n{pz) of the longitu- 
dinal momenta is then accounted for only in the end by 
averaging the results obtained with sharp values of Pz ■ 

A priori, such a treatment is no longer valid for a gen- 
eral grating interaction where different longitudinal mo- 
mentum components of the beam get correlated. We will 
accordingly use a two-dimensional scattering formulation 
to describe the grating interaction in Sect.lIVI However, 
we will see that for the parameters of typical experi- 
ments the main effect of taking the grating interaction 
beyond the eikonal approximation is on the transverse 
degrees of freedom. Since the effect on the longitudinal 
motion is much weaker we will retain the effectively one- 
dimensional description outlined above, postponing the 
physical discussion why this is permissible to Sect. lIVI 

Another idealization found in the basic treatments of 
the Talbot-Lau effect is to assume the transverse mo- 
tion in front of the interferometer to be in a completely 
incoherent state. This would correspond to a constant 
distribution of the transverse momenta p, and for gen- 
eral grating transformations, which depend on p, it is no 
longer a valid approximation. We will therefore present a 
formulation that takes into account a realistic beam pro- 
file and that allows for a general state transformation at 
the diffraction grating. Going beyond the idealization of 
a perfectly incoherent state will also permit us to assess 
the adjustment precisions required in an experimental 
implementation. 



A. The Wigner function and its transformations 

We start by briefly outlining how to describe matter 
wave interference by means of the phase space representa- 
tion of quantum mechanics |26l.l33l|. As discussed above, 
one may restrict the dynamics to the transverse beam 
state, which is most generally specified by its density ma- 
trix p. This state is equivalently described by the Wigner 
function 
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(1) 



which is a function of the transverse phase space coor- 
dinates {x,p). Like p it depends parametrically on the 
longitudinal momentum pz. Since we assume the beam 
to be coUimated w{x,p) is non-zero only for \p\ <C Pz- 



The great advantage of the phase space formulation is 
that it permits a straightforward, yet reahstic, descrip- 
tion of the beam and its propagation through the inter- 
ferometer. Most importantly, the free time evolution of 
a state during the time t is given by the same shearing 
transformation as in the case of the classical phase space 
density. 



wt (x,p) =wo\x p,p 



t 



(2) 



with m the particle mass. On the other hand, a quantum 
state transformation of the form p' = U/oU^, such as the 
effect of passing a grating, reads in phase space 

w'(x,p)= dxodpoK {x,p;xo,Po)w{xo,Po) ■ (3) 

Here the integral kernel is given by the the propagator 

K{x,p;xo,Po) = ^JJdsdsoe'^P'+P°'°^/^ 

x{x-^\U\xo+'^){x+'-\U\xo-'^y 



(4) 



Specifically, an ideal grating is characterized by a grat- 
ing transmission function t(x), with \t{x)\ < 1, de- 
scribing the multiplicative modification of an incoming 
plane wave. (The factor t{x) is non-zero only within the 
slit openings of the grating and it may there imprint a 
complex phase to account for the interaction potential 
between the grating walls and the beam particle, see 
Sect. IIIIl ') For such gratings the transformation reads 
(x|U|xo) = t (x) S {x — Xq), so that the grating propaga- 
tor reduces to a convolution kernel 



K{x,p;xo,Po) = 5{x-xo)^j dse'^P~P^>/'' 

x^(.-0r(. + |), (5) 

which is local in position [20|. This choice of K will be 
required below to reduce the generalized formulation of 
Talbot-Lau interference to the eikonal approximation. 

The convolution kernel provides a descriptive picture 
of the diffraction process. Suppose the incoming beam is 
perfectly coherent, i.e. a plane wave characterized by the 
longitudinal momentum p^ and transverse momentum po, 
corresponding to the (unnormalized) transverse Wigner 
function wq(x,p) — S {p — po). The diffracted transverse 
state is then given by 



Wl 
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(6) 

In case of a periodic grating with period d, the transmis- 
sion function can be decomposed into a Fourier series. 



iW= E 



bj exp 



{^nij 



(7) 



so that, after a free propagation over the longitudinal 
distance L, the transverse spatial density of the beam 
state is given by the marginal distribution 

P 



W2 (x) — I dpwi [ X L.p 



= '^ B,ni m— j exp \2Trim-] . (8) 
Here, we introduced the basic Talbot-Lau coefficients, 

oc 

BrniO^ E b,b*_^exp{nT^{m-2j)), (9) 



j = -oo 



and the characteristic length scale ix = (fi/X is called 
the Talbot length. 

Compare ^ to the Fourier coefficients of the trans- 
mission probability |i(x)p, which are given by the con- 
volution Am = J2i bjb'j-m- Equation ([5]) thus implies 
that the density distribution takes the form of the grat- 
ing transmission profile whenever the distance L is an 
integer multiple of the Talbot length, L — kL^, 



W2 (x) 



t[x + k- 



(10) 



j = -oo 



This is the elementary Talbot effect [3J], and it is the 
backbone of the Talbot-Lau interferometer, which how- 
ever does not require a coherent illumination. 



B. The general Talbot-Lau effect 

The general setup of a Talbot-Lau interferometer is 
described in Fig. [1] along with the most relevant param- 
eters. An incoherent particle beam represented by the 
shaded area passes a preliminary coUimation slit on the 
left side of the figure and enters the Talbot-Lau inter- 
ferometer through the first grating. The angular distri- 
bution of those particles in the beam which are finally 
detected is characterized by the spread a, determined by 
the collimation slit and the finite area of the detector. 

The first grating can be understood as an array of col- 
limation slits of period di preparing, after a distance L, 
the transverse coherence required for diffraction at the 
second grating. In this region, the distortion of a mat- 
ter wave front due to the grating interaction may signifi- 
cantly affect the interference pattern, so that the grating 
thickness b enters the calculation. (It is replaced by the 
laser waist Wz in case of a light grating.) 

The different diffraction orders still overlap further 
downstream if the distance is on the order of the Talbot 
length Lt — d"^ /\, and they may thus interfere among 
each other producing a near field fringe pattern. This 
can be explained qualitatively from the elementary Tal- 
bot effect ([To]) by considering each of the slits in the first 
grating as independent point sources. The Talbot pat- 
terns due to the individual slits will add up construc- 
tively for appropriate choices of the grating periods di. 
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FIG. 1: Schematic of a Talbot-Lau interferometer with the relevant setup parameters for a material diffraction grating. The 
beam emanating from the left is constrained by coUimation slits to an angular spread a before it enters the interferometer. The 
first grating of period di and size G S> di modulates the beam so that spatial coherence is created at the position of the second 
grating (z = 0) where the diffraction takes place. This leads to an interference pattern in the transverse beam density, ideally 
located at the position of the third grating. The latter can be shifted in the transverse direction by a variable distance xs- 
It thus modulates the flux provided its longitudinal position -qL and its period da match the interference pattern. A detector 
of size W further downstream measures the total incoming beam intensity as a function of xs- The central grating may be 
replaced by the phase grating created by a standing light field. 



d and of the distance factor 77, so that a distinct density 
pattern is created. It is verified without the need for a 
spatially resolving detector by superposing a third grat- 
ing whose period ds equals that of the density pattern, 
and by measuring the total flux through it as a function 
of its transverse position xs- 

We will now present a quantitative formulation of the 
successive propagation of the beam through the interfer- 
ometer. As discussed above, one can take the longitudi- 
nal motion of the beam particles to remain unaffected by 
the gratings. We may therefore assume the longitudinal 
momentum to have a definite value p^ for the time being, 
so that the longitudinal position z plays the role of a time 
coordinate t ~ zm/pz for the transverse motion. 



1. Sequential calculation 

The transverse state of the beam entering the interfer- 
ometer is far from pure. It is confined in position by the 
orifice of the first grating, whose size G typically covers 
thousands of grating slits. The momenta p are charac- 
terized by the angular distribution D (p/pz) whose char- 
acteristic spread is denoted as a. 

In a typical experimental situation, a rarely exceeds 
1 mrad corresponding to a fairly well-coUimated beam. 
However, this still covers a range of transverse momenta 
that is by orders of magnitude larger than the grating 
momentum h/di which separates the different diffraction 
orders. This explains why diffraction does not need to be 
taken into account at the first grating. The transverse 
Wigner function merely gets modulated by the grating 



profile, so that behind the first grating it reads 



wi {x,p) 



Gpz 



\ti{x)\'D[^ 



(11) 



Here ti (x) is the transmission function which is confined 
to the grating orifice, in principle. However, since the 
size G is typically larger than the sensitive region W of 
the final detector one may equally take it to be an un- 
constricted periodic function. 

The free propagation of the beam over the longitudinal 
distance L yields W2{x,p) = wi{x~ Lp/pz,p). The effect 
of passing the second grating is in general described by 
an operator U to be specified below. The corresponding 
phase space transformation is given by ^ so that, af- 
ter a second free propagation over the distance rjL, the 
transverse beam state in front of the third grating reads 

W3{x,p) = -— dxodpoD i—j \ti{xo)f 

xk(x- —vL,p;xo + ^L,p„) . (12) 

\ Pz Pz J 

The corresponding spatial density distribution, given by 
W3, {x) = J dpws {x,p), is now modulated by the third 
grating as a function of the grating shift 2:5 . The detec- 
tion signal is thus obtained as 



S{xs) 



W/2 



dxws (x) 1*3 (x- Xs)\'' 



w/2 

CXD 

2J Sn exp f 27rm— J , 



(13) 



where jts (a;)| is the spatial transmission probabihty of 
the grating and W the size of the detector. The latter 
typically covers hundreds of grating periods so that one 
may disregard the finiteness of W when evaluating the 
basic interference effect. However, as shown below, it 
does play a role when the experimental adjustment pre- 
cision needs to be evaluated. 

The resulting interference pattern is most easily char- 
acterized by the contrast of the modulation signal, con- 
ventionally defined as the ratio between the amplitude of 
the signal variation and its mean value. However, since 
the experimental signal is usually very close to a sinu- 
soidal curve, it is in practice more convenient and more 
robust to use the sinusoidal visibility defined in terms of 
the first to the zeroth Fourier expansion coefficient of the 
da-periodic signal [25| . 



V = 



2Si 



(14) 



It can be easily obtained from noisy experimental data by 
fitting a sine curve and it coincides with the conventional 
definition in case of a sinusoidal signal. 



2. Decomposing the grating propagator 



which are the central quantities for describing the inter- 
ference effect. 
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X exp (iTT^ (m — 2j)) , 



(18) 



We note that at integer values of the argument ^ the 
phase factor in (|18p reduces to a constant sign, while p 
has the meaning of an incident momentum and i> that of 
a scale factor, as will be seen below. 

By using the generalized Talbot-Lau coefficients the 
density distribution at the third grating takes the explicit 
form 



wsix) 



1 E A,,e^J2m(rn + n^]^ 



m,.n—~oo 



dp f pL ( d 

X / — exp —ZTTi rim -I- (f; -I- l)n — 

Pz \ Pzd V "1 

(19) 

This expression for the interference pattern is now fur- 
ther simplified by identifying the Talbot-Lau resonance 
condition. 



We proceed to evaluate the interference pattern ([T^ by 
noting a general property of the grating transformation 
U. The periodicity of the grating implies that the po- 
sition representation (a;|U|a;o) is d-periodic with respect 
to the center position (a; -I- a;o)/2. This admits the series 
expansion 



(x|U|xo) 



E ^^P 



mn ; — 1 Un{x — xa) (15) 



where the transformation within a single slit is now char- 
acterized by the corresponding Fourier coefficients 



Un {x) = - 
d 



d/2 



d/2 



1- f r. ^ 

dx exp —ZTTtn— 
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f|U|. 



x^ 
2' 



(16) 

It is now convenient to specify these transformation func- 
tions in terms of their shifted Fourier transformation, 
which serves to generalize the grating coefficients intro- 
duced for the eikonal case in ^. 



I' 



nh 



bn (p) ^ I dx exp { -T \P+ ^r" ] x] ^n (x) (17) 



Indeed, for position-diagonal operators U, these functions 
drop their momentum-dependence and reduce to the co- 
efficients in d?]). 

The generalized grating coefficients (1171) can now be 
used to construct the generalized Talbot-Lau coefficients. 



3. The resonance condition 

Only a part of the double summation in ()19p con- 
tributes appreciably to the interference pattern. This 
follows from the fact that the ratio L/d of grating dis- 
tance to grating period is a large number, typically on 
the order of 10^. As a result, the momentum dependence 
of the phase factor in P^ occurs on a very different scale 
compared to the variation of the momentum distribution 
D{p/pz) and of the Talbot-Lau coefficients Bm{^;p,i^)- 
In fact, in the idealized case of both a completely incoher- 
ent illumination and an eikonal interaction the latter two 
functions are independent of p, so that the momentum 
integral is finite only if the phase vanishes identically, i.e., 
for 



rjm 



{r] + l)—n 
di 



0. 



(20) 



It would imply that an interference pattern is observed 
only if (?7 + \)d/{rjdi) is a rational number. 

This strong resonance condition gets relaxed if we ac- 
count for the weak momentum dependence of the remain- 
ing integrand in ([TO]) . We assume that those index pairs 
(to, n) of the double sum contribute appreciably where 
the phase variation in the momentum integration does 
not exceed about tt. Since the value of \p\/Pz is bounded 
by the angular spread a this leads to the relaxed condi- 
tion 



V 



1 d 



7] di 



J_d 
- 2^L' 



(21) 



Here r and s are natural numbers without common divi- 
sor. They indicate the type of resonance and specify the 
set of index pairs (m, n) e {{r£, si) : £ Q Z} contributing 
to the sum. 

If a is about Imrad and L/d = 10'°^ the right hand 
side of (PT|) is on the order of 10^^, and in practice only 
a single resonance dominates for each set of parameters 
d, di, r]. The expression for the interference pattern thus 
simplifies to 
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Here, the period of the interference pattern is given by 



*^l5-il ■ <^'' 



Note that a large interference contrast is obtained for the 
low order resonances (with small values of the integers r 
and s) because the Fourier and Talbot-Lau coefficients 
Arm Bm gcnerically decrease with their order m. 

The standard choice in experiments is to take equal 
gratings, d = di, in an equidistant configuration, -q — 1. 
This corresponds to the r : s = 2 : 1 resonance. One 
may obtain even larger contrasts with 77 = 1 by using the 
basic r : s = 1 : 1 resonance, at the expense of dealing 
with different gratings, di — 2d. 

Also magnifying or demagnifying interferometers can 
be realized [2^]. Consider, for example, the case where 
the first two gratings are equal, d = di. The period of 
the interference signal is then given by d^ — d/(r — s) 
and cannot thus be greater than d. The condition 
r/s = {r] + 1) /rj with 77 > implies that r must be 
strictly greater than s. If one wishes to decrease da, one 
must choose r — s > 2 which requires rj < 1, i.e., a com- 
pressed setup. However, the values of r and s grow with 
decreasing ^3 so that the visibility of the resulting in- 
terference pattern decreases significantly. A magnifying 
interferometer, on the other hand, can only be realized 
with different gratings d ^ di. 

If the aim is to create a high contrast interference 
pattern with a specific period ^3 a low-order resonance 
should be taken, such asr:s = 2:lorr:s = l:l. 
The experimental setup parameters {ri,d,di) must then 
satisfy the equations 



{7] + l)d 
rjdi 



rds 
sds 



(24) 
(25) 



with desired resonance parameters (r, 5,^3). The solu- 
tion is not unique, in general, so that one has a certain 
freedom to account for experimental limitations. 



C. The Talbot-Lau effect in eikonal approximation 

The expression for the interference pattern can be fur- 
ther simplified if the grating interaction is treated in 
eikonal approximation. The grating coefficients (|17p turn 
then into the momentum-independent Fourier coefficients 
bj of the transmission function ([7]). Similarly, the gen- 
eralized Talbot-Lau coefficients (TTSt then reduce to the 
basic coefficients given in Eq. ^ so that the prediction 
for the density pattern P^ simplifies to the series 
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(26) 



It involves the Fourier transformation of the angular dis- 
tribution, 



J Pz \PzJ 



(27) 



In the case of a very broad momentum distribution this 
characteristic function D can be replaced by a Kronecker- 
5 function, which yields the basic result [2^, [26[ 
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The comparison with Eq. ([8]) shows clearly that Talbot- 
Lau interference is based on the Talbot effect, both de- 
scribed by the Talbot-Lau coefficients defined in ^. 

It should be emphasized, though, that in order to es- 
tablish that an observed fringe pattern is really due to 
a quantum effect one must compare the quantum pre- 
diction with the result of the corresponding classical cal- 
culation. This is necessary since a classical moire-type 
shadow pattern might give rise to a similar observa- 
tion, even tough the classical contrast is typically much 
smaller. It is shown in the appendix how the classical 
treatment can be formulated in the same framework as 
the quantum case using a phase-space description. 

Finally, since Eq. (P5|) assumes the resonance condition 
to be exactly met, it cannot be used to assess the adjust- 
ment precision required in the experiment. As shown 
in the next section, one can evaluate the necessary pre- 
cision in the eikonal approximation by explicitly taking 
into account the finite transverse momentum spread and 
the finite size of the signal detector. 



D. Adjustment requirements 

In order to allow for deviations from the exact reso- 
nance condition we use Eq. (|26p when evaluating the ex- 
pression for the detection signal. According to Eq. (fT^ 
the detection signal S{xs) is then characterized by the 
Fourier coefficients 



t>n 



w 
"g 



A': 



E 



Af B„ 



di 



L 



xD 



2t:L ( —m 
a 



V 



1 



X suic 



a 



(29) 



di 

£ n 
di rfa. 

Here, the A'^ are the Fourier expansion coefficients of the 
third grating proffie jia (x) \ with the period ^3. Consider 
now small deviations {SL, Sr], 5di,Sd, Sds) from the setup 
parameters (L, 77, di, d, ds) satisfying a particular Talbot- 
Lau resonance condition r : s. The distances Li = L and 
L2 ~ 1]L can thus vary independently. Instead of approx- 
imating the characteristic function Z? by a Kronecker-5, 
we account for the small deviations by the refined ap- 
proximation D {uj + e) w S^flD (e). It holds as long as 
the main argument uj is either zero or much greater in 
modulus than the width l/a, and as long as e < \/a. 

One can thus split the argument of D in (j29|l into the 
ideal resonance part lo and in the part e containing the 
parameter deviations. The same procedure may be done 
with the sine term, which is sharply peaked if the detector 
size W exceeds the period d by orders of magnitude. If 
the small parameter deviations are taken into account to 
first order one obtains the same resonance relation for 
the summation indices in (|29p as in the ideal case, while 
the nth Fourier expansion coefficient of the signal (with 
respect to the period d^) is multiplied by the reduction 
factor 



Rn 



D 



( l-KsnL 



\ di 



5L_ 
T 



-■^ + iv + i) 



X smc 



TTsnW 
di 



Srj 

V 
Sdi 

~d^ 



Sd 
~d 



5d 
~d 
6d^ 
ds 



6(h 
~d^ 
Sd' 
~d 



(30) 

The absolute value of this factor is less than 1 for n ^ 
so that the contrast of the detection signal is effectively 
reduced by the deviations. In particular, the sinusoidal 
visibility P^ of the signal is reduced by the factor Ri. 
Typical experiments [20, 1211, 122, l25l are characterized by 
W/d « 10^, L/d « 10^ and a fv lO'^ so that aheady rel- 
ative parameter deviations on the order of 0.1% strongly 
affect the interference contrast. 

To obtain a simple and conservative estimate we treat 
all imprecisions as independent, specializing to the stan- 
dard case r : s = 2 : 1, where the grating distances 
and grating periods are equal. Bounds for the adjust- 
ment precisions are then obtained from Eq. pop by re- 
quiring the arguments to be smaller than the widths of 



the momentum distribution and of the sine function, re- 
spectively. 



and 



SL 
L 


+ 2 


Sd 
~d 


d 
AnaL 




Sd 




d 




~d 


<, 


AW' 



(31) 



(32) 



This quantifies to what degree a better collimation of 
the beam and a smaller detector size relax the required 
adjustment precision, albeit at the expense of a loss of 
signal. 



III. EFFECT OF THE INTERACTION 
POTENTIAL 

The preceding section showed how the general coher- 
ent state transformation effected by a diffraction grat- 
ing enters the Talbot-Lau calculation via the generalized 
grating coefficients (IT71) . In general, this transformation 
is determined by the potential V (x, z) due to the long- 
range dispersion forces acting on the beam particles while 
they pass the grating structure. 

Before we present a general way to account for the 
presence of this potential in Sect. IIVI it is helpful to dis- 
cuss the most important grating-particle interactions in 
a simpler form, by using the eikonal approximation. It 
treats the grating as the combination of an absorption 
mask and a phase modification, and it can be character- 
ized by a grating transmission function 



t (x) = \t {x)\ exp I -■— dzV (x, z) 1 , (33) 

whose amplitude |t(-)| '-^ {0,1} describes the grating 
structure. Interaction-free gratings are modeled by a 
grating transmission function without phase, t{x) — 
|i(x)|, while pure phase gratings, such as the standing 
laser field, wave are characterized by \t{x)\ = 1. The cor- 
responding eikonal interference pattern is then obtained 
immediately as described in Sect. Ill CI 

We proceed with a short overview of the typical grat- 
ing potentials, discussing how they affect the interference 
contrast in the eikonal approximation. By convention, 
the diffi-action grating is located at the longitudinal po- 
sition z = 0, as indicated in Fig. [T] 



A. Material gratings 

A neutral particle located within the slit of a mate- 
rial grating experiences a potential determined mainly 
by the attractive dispersion forces due to the grating 
walls. Other forces, such as the exchange interaction or 
the electro-static attraction due to a permanent dipole, 
are much less important for interferometry because they 
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0.45 




4 5 6 7 8 9 

de Broglie wave length A, (pm) 

FIG. 2: Talbot-Lau visibilities ([14} for the equidistant setup 
with three material gratings [20| as a function of the de Broglie 
wavelength of the interfering C70 fuUerenes. At the second 
grating we assume no dispersive interaction (dotted line), a 
van der Waals interaction (dashed line), and a retarded van 
der Waals interaction (solid line) . The dash-dotted line repre- 
sents the classical calculation with a retarded interaction. It 
depends on the fictitious "wavelength" A = h/{mvz) via the 
velocity v^ because the interaction time is velocity dependent. 



either occur only at very close distances or because they 
are diminished by rotational averaging. In any case, for 
all positions within a slit the walls are well approximated 
by an infinite surface in the yz-plane [33,[3a| • The grating 
potential is thus set to be solely x-dependent and acting 
only within the time of passage t = mb/pz throughout 
the grating of thickness &, while it is set to be zero out- 
side. 

In general, the dispersion force between a polarizable 
particle and a material plane is described by the expres- 
sion of Casimir and Polder [2^ and its generalizations, 
e.g. [33, [33. In the close distance hmit it reduces to 
the van der Waals potential V{x) = —C^/x^, with x 
the distance to the surface, and at large distances, where 
retardation plays a role, one has the asymptotic form 
V{x) = —C^jx"^. The interaction constants Cj,^C^ > 
are determined by the frequency dependent polarizabil- 
ity of the particle (as well as the dielectric function of the 
grating material) , and the regime of validity of the limit- 
ing forms is delimited by the wavelengths corresponding 
to the strong electronic transitions [Sg, [33, EQ] . 

In any case, the dispersion force diverges on the grat- 
ing walls, rendering the eikonal approximation invalid in 
close vicinity to the walls. Since the contributions of 
these regions do not alter the interference contrast ap- 
preciably, but lead to numerical noise, we will discuss a 
reasonable criterion to blind out the beam close to the 
grating walls in Section ITVl 



Figure [2] shows the interference visibility (fT4)) for a typ- 
ical Talbot-Lau experiment with C70 fuUerene molecules 
of mass m — 840 amu [20] , plotted versus their de Broglie 
wavelength. All gratings are assumed to be separated by 
the distance L = 22 cm and to be made of gold with a 
period oi d — 991 nm, a slit width of 476 nm, and a thick- 
ness oib = 500 nm. One can see two peaks corresponding 
to the first and second Talbot order, at A = 4.5 pm and 
at A = 9.0 pm, respectively. The solid line represents 
the eikonal approximation with a retarded asymptotic 
potential (C4 — Shcao/Sir with the static polarizability 
ao = 96.7 A^ obtained via the Clausius-Mossotti rela- 
tion [41|), while the dashed and the dotted lines corre- 
spond to the van der Waals potential (C3 = 10 meVnm^) 
and to the absence of an intra-slit potential, respectively. 
One observes that the presence of the dispersion forces 
changes the interference characteristics significantly. The 
asymmetry in the double-peak structure compared to the 
interaction-free case is due to the fact that the particles 
with a larger velocity, i.e., with a smaller wavelength A, 
receive a smaller eikonal phase than the slower particles. 
Moreover, note that the effect of retardation has a small 
but visible influence on the fringe contrast. 

The dash-dotted line in Fig.[2|represents the moire-like 
effect as expected from classical mechanics. It was calcu- 
lated in a phase space formulation as described in the ap- 
pendix, using the classical correspondence of the eikonal 
approximation with the C4-potcntial. As one expects, 
the dependence of the classical result on the fictitious 
"wavelength" A = h/{mvz), which is due to the veloc- 
ity dependence of the classical deflection, differs strongly 
from the quantum results. 

We emphasize that all numerical results presented in 
this paper are obtained for a fixed longitudinal velocity 
Vz = h/{Xm) of the particles. A comparison with the 
experiment still requires the results to be averaged with 
respect to the velocity distribution in the beam. This 
may pose a severe restriction when using particles with 
larger polarizability because the increased dispersive in- 
teraction decreases the width of the double-peaks seen in 
Fig. [H significantly ^^ . One way to avoid this is to re- 
place the material diffraction grating by a standing laser 
wave. 



B. Laser gratings 

Recently, a Kapitza-Dirac- Talbot-Lau interferometer 
(KDTLI) was demonstrated, where the central grating 
is formed by a standing light wave [22]. In the eikonal 
approximation the phase of an incoming plane wave is 
modulated according to the potential created by the off- 
resonant interaction with the standing laser beam. It is 
determined by the energy of the induced electric dipole 
in the oscillating field, and is therefore proportional to 
the laser power Pl and to the dynamic polarizability a^ 




de 



1.5 2 2.5 3 

Broglie wave length X (pm) 



FIG. 3: Fringe visibilities (fT4)l for a Kapitza-Dirac-Talbot- 
Lau interferometer [22| in the eikonal approximation. The 
solid line gives the quantum result for C70 fuUerenes, while the 
dashed line represents the corresponding classical calculation 
(neglecting photon absorption). Here we assume a laser power 
of Pl = 6 W, a vertical waist of Wy = 900 fim, a grating 
distance of L = 105 mm, a grating period of d = 266 nm, and 
a dynamic polarizability a^ — 118 A'^ corresponding to a laser 
wavelength of 532 nm. 



([35)) . In contrast to Figure [2l the visibility drops to zero 
whenever the wavelength corresponds to an integer multi- 
ple of the Talbot condition. This is explained by the fact 
that in the elementary Talbot effect PU)) a pure phase 
grating with no absorptive walls leads to a constant den- 
sity at multiples of the Talbot length. At the same time, 
there are broad regions of high contrast. They render the 
KDTLI setup superior to a material grating interferome- 
ter for particles with a high polarizability and a substan- 
tial velocity spread. Moreover, the interaction strength 
and the grating period can be tuned in a KDTLI render- 
ing the laser grating the preferable choice to explore the 
validity regime of the eikonal approximation, as done in 
Section El 

The dashed line in Figure [3] represents the classical 
version of the eikonal calculation as described in the ap- 
pendix. Note that it differs significantly from the quan- 
tum result in the experimentally relevant wavelength 
regime, while the curves become indistinguishable for 
A — > 0. Finally, it should be emphasized that a real- 
istic description of the laser grating must also account 
for the possibility of photon absorption. The effect of 
the resulting transverse momentum kicks can be incor- 
porated into the eikonal approximation by replacing the 
grating transformation kernel by a probabilistic sum of 
such kernels 1231. 



of the beam particles at the laser frequency. 



V{x,z) 



4:PLauj 
TreQCWyWz 



^^"' ("?) 



-2zV»' 



(34) 



Here, 



and Wz are the waists of the Gaussian mode. 



and Wy is chosen large compared to the detector size in 
order to guarantee a regular grating structure. We there- 
fore disregard the y-dependence by setting y = 0. More- 
over, for sufhciently small Wz an effective one-dimensional 
treatment of the potential is permissible, where one re- 
places the z-dependence by a parametric time depen- 
dence z — pzt/m, as discussed below in Sect. IIVCI 

The z-integration in (j33p renders the eikonal phase in- 
dependent of Wz, which already indicates that the ele- 
mentary eikonal approximation will cease to be valid if 
the laser waist is increased. Nonetheless, it yields an in- 
terference contrast that fits well to the measured data of 
the recent fullerene experiments |22l |. Moreover, it admits 
a simple, closed expression for the Talbot-Lau coefficients 



Bm (C) — Jrn —- 



'QjrhceQWyP; 



■ sin 7rf 



(35) 



where M is the mass of the beam particles and Jm stands 
for the mth order Bessel function of the first kind 42] . 

The solid line in Figure [3] shows the visibility P^ for 
the Viennese KDTLI with fuUerenes 221 as obtained from 



IV. SEMICLASSICAL APPROACH TO THE 
GRATING INTERACTION 



The de Broglie wavelength is by far the smallest length 
scale in the matter wave interference experiments consid- 
ered in this paper. It seems therefore natural to use a 
semiclassical approach for calculating the grating trans- 
formation used in p6)) . We will show how the elementary 
eikonal approximation (j33p can be derived from the ap- 
propriate semiclassical formulation if a high-energy limit 
is taken. Moreover, one can obtain a refined eikonal ap- 
proximation, where an incoming plane wave is still merely 
multiplied by a factor. Since this factor depends on the 
transverse momentum of the incoming wave it renders 
the numerical implementation more elaborate. 

However, the main result of this section is the expres- 
sion ([M)) - ([55)) for the scattering factor which approxi- 
mates the scattering transformation of a transverse plane 
wave in the semiclassical high-energy regime. As such it 
permits to evaluate the generalized grating coefficients 
P?)l straightforwardly, see Eq. ([55)) . As shown in Sec- 
tion |Vl this result outperforms the eikonal approxima- 
tion and it may serve to characterize its regime of va- 
lidity. The derivation is based on the semiclassical ap- 
proximation of the two-dimensional time evolution oper- 
ator Ut determined by the grating interaction potential 
V{r) = Vix,z). 
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A. The scattering factor 



B. Semiclassical calculation 



We proceed to incorporate the grating interaction by 
means of the formahsm of scattering theory [43| . Its basic 
tool is the scattering operator defined as 



S= hm UL",^U2*UL°], 

t—>-oo 



(36) 



with Uj = exp {^—ip^t/2mfij the free time evolution op- 
erator for the motion in the xz-plane and Uf the com- 
plete time evolution operator, which includes the grat- 
ing potential V{r). It is pertinent to use S instead of 
lit since the scattering operator transforms the state in- 
stantaneously leaving the asymptotic dynamics to be de- 
scribed by the free time evolution. Since the latter is 
easily incorporated using Wigncr functions this fits to 
the phase space description of Sect. IIIBl where the ini- 
tial beam state entering a Talbot-Lau interferometer is 
propagated freely to the second grating before the scat- 
tering transformation is applied. In the expression ^ 
for the grating propagator, which serves to calculate the 
interference pattern (|22p , one may thus use the S-matrix 
([55)1 in place of the general unitary operator U. 

In the basis of the improper plane wave states {r\p) = 
{2TTh)~^ exp(ir • p/h) the scattering operator ([36]) is con- 
veniently described by the scattering factor 



$(r,p) 



(r-|S|p) 
(r\p) 



(37) 



Notice that S acts in the Hilbert space defined on the 
two-dimensional plane r — {x,z). However, as will be 
justified below, the longitudinal motion may be separated 
and treated classically so that the transformation is con- 
fined to the transverse dimension, while the z-coordinate 
turns into an effective time coordinate for a given lon- 
gitudinal momentum p^. The brings about the reduced 
scattering factor 



(j){x,p) = $(xea;,pe^ +Pzez) 



(38) 



which depends parametrically on pz and which is evalu- 
ated at the position z = of the center of the diffraction 
grating. 

The reduced scattering factor (j>{x,p) describes the 
phase and amplitude modification of a transverse plane 
wave with momentum p due to the grating interaction. It 
enters the Talbot-Lau calculation via the Fourier coeffi- 
cients (J16p after an expansion in the plane wave basis. It 
follows that the generalized grating coefficients (fT7|) are 
directly related to {x,p) by a Fourier transformation. 



bn (p) 



d/2 



d/2 



da;e-2""=^/''<?!)(a:,p). 



(39) 



The calculation of the Talbot-Lau interference thus re- 
duces to evaluating the scattering factor 



We proceed to calculate the scattering factor (pS]) by 
means of the semiclassical asymptotic approximation. It 
assumes the action of those trajectories through the in- 
teraction region, which contribute to the path integral 
for its position representation, to be much larger than 
h, and the particle wavelength to be much smaller than 
the scale where the interaction potential changes appre- 
ciably, |AVy| ^ V). These conditions are well satisfied 
in the typical experimental situation if we disregard the 
regions very close to the gratings where the potential ex- 
ceeds the kinetic energy. We may thus approximate the 
time evolution operator in the position representation by 
the semiclassical van Vleck-Gutz wilier propagator [4J] 



(r-|Ut|ro) = 



1 



2nih 



det 



/ d^St{r,ro) 
\ drdro 



X exp ( -S't (r,ro) 



(40) 



Here, St is the action of the classical trajectory travelling 
during time t from the position Vq to r. In general, there 
might be more than one such trajectory, which would re- 
quire taking special care of almost coalescing trajectories 
and of the associated Morse index. However, in the inter- 
ferometric setup we are in the high-energy regime, where 
the interaction potential is much weaker than the energy 
of the incoming particles, \V\ <C E. All relevant contri- 
butions are therefore characterized by a single, slightly 
deflected classical trajectory passing the interaction re- 
gion. 

It is now convenient to specify this trajectory in terms 
of the deviation from the undeflected straight line, as 
specified by the initial position Tq and momentum Pq. 
The momentum change after time t is given by 

Ap,(ro,Po) = - / dTVF(r,(ro,Po)), (41) 
Jo 

so that the deflected trajectory reads 

rt{ro,Pa) = To + — Po + / — Ap^(ro,Po) (42) 



Pt{ro,Pa) = PQ + '^Pt{ro,Po) 



(43) 



In the van Vleck-Gutzwiller propagator (|^(I)) the con- 
tributing trajectory is specified by the boundary values 
ro and r. For the following calculation it is important 
to rewrite it as an initial value problem, specified by the 
initial phase space point (ro,Po) of f^^ trajectory |4a |. 



(^|Ut|ro) 



1 



2nih 



d Po- 



det 



/ 9rt(ro,Po) 
V 5Po 



xJ(r-rt(ro,Po))exp ( -St{ro,Po) 



(44) 
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Plugging this into the expression ((36| for the 2d scatter- 
ing operator yields a semiclassical approximation for the 
scattering factor ([37| . 



$(r,p) 



rn 
: lim — - 

T^oo T 






X exp T©2T(ro,Po) 



det 



/ 9r2T(T'0,Po) 

V 5Po 



(45) 



The phase of the integrand is given by the action- valued 
function 

02T(r-o,Po) = '5'2T(ro,Po)- ^(^-^2T(r-o,Po))^ 

(46) 



-^•("-""-£^^, 



while the amplitude is determined by the stability deter- 
minant of the associated trajectory. 

The semiclassical van Vleck-Gutz wilier propagator 
()40p used here is a stationary phase approximation of 
the time evolution operator in path integral formulation 
[4J] . To remain at a consistent level of approximation it 
is therefore necessary to evaluate the phase space integral 
in (P5|) in the stationary phase approximation as well [4g| . 
The stationary point of the phase 82T is determined by 
the condition 



de 



2T 



dro 



de 



2T 



dPo 



0. 



(47) 



Noting the initial value derivatives of the classical action 

T / \ T 





Pt 



(48) 



(49) 



and using the fact that the matrix dr2T/dpQ is invertible 
for our trajectories, the stationary phase condition leads 
to the equations 



r 
P 



r2T{ro,Po) 
Po- 



T 



P2Tiro,Po) 



(50) 
(51) 



They serve to determine the initial position ro{r,p) of 
the trajectory implicitly. Using the general formula for 
a four-dimensional integral, e.g. [43, Eq. (A. 30)], the 2d 
scattering factor (|45p then takes the form 



'^ir,p) = lim A2T(ro(r,p),p)exp -Q2T {roir,p) ,p) 

(52) 
The amplitude modification 

1 
'dr2T{ro,Po) m dp^T iro,Po) 



^2T(r-o,Po) 



dot 



dro 



T 



dro 




(53) 



t=2T 



FIG. 4: Sketch of the stationary phase conditions for the semi- 
classical scattering factor $(r,p) defined in (I52|) . While the 
momentum p fixes the asymptotic initial momentum p^ of the 
classical trajectory, the initial position vq is determined im- 
plicitly by the spatial point r. The latter defines the asymp- 
totic final position r2T of the trajectory by means of a free 
evolution determined by the asymptotic final momentum P2T- 
For interaction potentials that admit a scattering theory de- 
scription this construction guarantees that the action (|46[) and 
the stability amplitude (I53|l converge as T — > 00. The range 
of the scattering potential is indicated by the shaded area. 



is obtained, in a tedious but straightforward calculation, 
by using the Poisson relation between the conjugate vari- 
ables Vt and Pf. of the trajectory [4g| when evaluating 
the product of two determinants. The matrix of deriva- 
tives in (j53p can be computed by taking the initial value 
derivative of the equation of motion for the trajectory 
and solving the resulting ordinary differential equation. 

Figure [4] shows how the stationary point can be un- 
derstood from the point of view of a classical scattering 
trajectory. Given the phase space coordinates r,p deter- 
mining the matrix element of the scattered plane wave, 
the momentum p fixes the initial momentum of the clas- 
sical trajectory passing through the interaction region. 
The position r determines the final position r2T of the 
deflected trajectory, which is obtained after a free motion 
during time T in a direction given by the final momen- 
tum P2T of the trajectory. This free evolution ensures 
that the associated action (|15)) and stability amplitude 
((53)) is independent of T in the limit T — > 00. Note that 
Fig. 0] overemphasizes the effect of the potential, since in 
all interferometrically relevant situations the deflection of 
the classical trajectories will be only a small correction 
to the free rectilinear path. 

The limit in the semiclassical result (|52p is obtained 
already at a finite time T if the interaction potential may 
be considered to have a finite range. The relevant scale is 
the grating passage time, given by i = mb/pz for material 
gratings and t = mwz/pz for laser gratings, respectively. 
This follows from the fact that outside of the grating 
potential the classical trajectories remain rectilinear, so 
that neither the stationary phase condition (|50p nor the 
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stationary phase 02t and amplitude A2T are affected by 
a further increase in T. 



a quite different argumentation [28|. The reduced scat- 
tering factor ([55]) for the transverse dimension then reads 



Iterative solution in the momentum deflection 



(t>ix,p) =exp 



-- — / azV \ X ~\ z, z 

n-Pz J -00 \ Pz 



(57) 



We wiU now give expUcit approximate solutions for the 
semiclassical scattering factor (|52p which are valid in the 
high-energy limit, i.e., whenever the longitudinal kinetic 
energy pl/2'm of the incoming particles is large compared 
to the interaction potential |y(r)|. In general, one has 
to solve Eq. (|50p for the initial value tq- Rewriting it in 
terms of the expressions (H^ and ([^5]) for the positions 
and momenta of the trajectory, one obtains an implicit 
equation for ro(T",p) which may be solved iteratively. 



T T 

ro{r,p) =r p + —Ap2T {ro {r,p) ,p) 

m m 

r2T 



Jo m 



The Glauber eikonal approximation 



(54) 



In the semiclassical short interaction time and high- 
energy limit, which applies to typical interference exper- 
iments, the transverse momentum deflection Ap2T i^ so 
weak that its contribution can be neglected. This corre- 
sponds to the zeroth order solution of Eq. ([54]) , 



,(0) 



T 

(r,p) =r p. 

m 



(55) 



In order to keep the approximation consistent one has 
to neglect the deflection in the expressions for the phase 
modification P5|) and for the amplitude modification ([55]) 
as well, by setting the classical trajectory to be the free 
rectilinear path. Consequently, there is no amplitude 
modification and the scattering factor ([5^ reads 



(56) 



$ (r,p) = exp — 



dtV {r + —t 
m 



This is the Glauber eikonal approximation obtained by 
R. J. Glauber from the Lippmann-Schwinger equation in 



The elementary eikonal approximation (|33|) follows from 
this expression in the limit of vanishing transverse mo- 
mentum, p — > 0. It applies in the case of a well-coUimated 
beam and a small grating thickness b, i.e., \pb/pz\ ^ d. 



2. The deflection approximation 

The Glauber eikonal approximation ceases to be valid 
as the interaction strength or interaction time increases, 
and one has to go to the first order in the momentum de- 
flection Ap( when evaluating the stationary initial value 

dMD, 



.(1) 



T 



T 



{r,p) = r pH Ap 



m 

r2T 



2T 



T 

—p,p 

m 



f di , / T 

/ — Apf r p, p 

Jo m \ rn 



(58) 



The higher order terms, which are neglected here, involve 
derivatives of the momentum deflection ApjT- 

The time-evolved trajectory starting from this initial 
value is approximated, again to first order in Ap^, by 



'^t {r''„^\'r,p),p 



T 



-P,P 



T 



Ap 



2T 



T 
r P,P 



(59) 



2T 



— Ap^ r p,p 

Q m \ m 



This trajectory must be used when calculating the action 
(|^5)l in order to ensure that all expressions are evaluated 
to the same order in the deflection Ap^ . Since this holds 
also for the time integral J dtV {rt (r"o,p)) a 1st order 
Taylor expansion of the potential is required. In total 
this yields 



J 



B2T(r-,p) 



2T 



dtV {rt{r 



-Ap. 



2T 



T 



-P,P 



T 



-P,P 



2T 



1 

2m ,/Q 



2T 



dtApl 



— Apt [r p,p 



T 



-P,P 



T 
2rn 



Ap 



T 
r P,P 



(60) 



On the other hand, when evaluating the amplitude (|53p the zeroth order solution of the initial value (j55[) must be used 
instead of ([55)1 because the classical equation of motion for the matrix of initial value derivatives {drt/drQ, dp^/dro) 
is governed by the derivative of the interaction force rather than by the force itself. Accounting for the modification 
in calculating the stability determinant would therefore amount to a higher order correction in the momentum 



13 



deflection. It follows that the amplitude factor is consistently approximated by 

/ dr2T {i — pT/m, p) T dp2T (^ — pTjm, p) 



A2T{r,p) 



det 






9ro 



TO 



Otq 



-1/2 



(61) 



r 



3. Separation of the longitudinal motion 

Based on the deflection approximation of the station- 
ary scattering factor we may now derive a reduced scat- 
tering factor ([38]) that can be incorporated into the one- 
dimensional Talbot-Lau calculation and that significantly 
extends the validity regime of the eikonal approximation. 
It is necessary for this purpose that the longitudinal mo- 
tion remains unaltered. To a very good approximation 
this is indeed the case for the small trajectory deflections 
required above in the deflection approximation. This is 
due to the energy conservation during the scattering pro- 
cess, ensuring (p -|- Ap2t)^ = P^ where p = {p,Pz) is the 
incoming momentum and Ap2T the total deflection (|4ip . 
It follows that a small transverse deflection Ap2T <C Pz 
yields in a well-collimated beam, \p\ ^ p^, a total lon- 
gitudinal deflection Apz,2T ~ p/PzAp2T which is much 
smaller than Ap2T- The longitudinal part of the classical 
trajectory ([M]) may thus be treated as a free motion with 
constant momentum pz , 



zt 



(roHr,p) 



,+ Pl^t-T) + o(Ap2T- 
m \ Pz 



This reduces all the vectorial quantities in the semiclas- 
sical phase and amplitude modification to the transverse 
scalar quantities. 

In addition, one can now remove the longitudinal mo- 



J 



tion altogether from the scattering description by switch- 
ing into a comoving frame with velocity Vz — Pz/rn. It is 
convenient to redefine the transverse components of the 
classical transverse trajectory so that they start at — T, 



^Pt{x,p) 



At 



-T 



d^v{xr{x,p),^T] (63) 



Xt (x,p) — X ^ 1 



Pt{x,p) = p + Apt{x,p). 



dr 



Apr{x,p) (64) 
(65) 



These are the scalar analogues of (|^ - (P5)) but for the 
shift in the time coordinate by —T, which accounts for 
the asymptotic initial condition of the free longitudinal 
trajectory ([S^ . 

We can now state the result for the reduced scattering 
factor (f>{x,p) = $ {xex,pex + Pz^z) to first order in the 
momentum deflection. 



(66) 



0(x,p) =a{x,p)exp ( -9{x,p) 



The amplitude is given by 

a{x,p) 
while the phase can be simplified as 



lim 



T 

dxXT {x,p) dxPT {x,p) 

m 



(67) 



0{x,p) 



lim 



At 



-T 



V ' 



{xt{x,p),—t] dxV {xt{x,p),—t] / ArdxV {xr{x,p),—T] 

\ m / m \ m / , \ my 



(68) 



r 



Here xt{x,p) is the classical one-dimensional trajectory 
of a particle starting, at t = — T, with momentum p 
at the position x — pT/m. It evolves in the effectively 
time-dependent potential Ut{x) — V {x,tpz/m). The as- 
sociated momentum is pt{x,p). It is easy to see that the 
limits in ([S7)) . ([55)1 exist for sufficiently short-ranged po- 
tentials. In particular, if the scattering potential has a 
finite extension the limit is reached already at finite times 
T, once the longitudinal distance 2Tvz is larger than the 
size of the interaction region. 

The reduced scattering factor ([66]) constitutes a signif- 
icant improvement over the eikonal approximation used 
so far, as will be demonstrated in the next section. Com- 
pared to the elementary eikonal approximation (|33p and 



to the Glauber eikonal result (l57|) . this expression not 
only provides the consistent incorporation of the deflec- 
tion into the phase, but it also introduces an amplitude 
modification of the scattered wave. 



V. NUMERICAL ANALYSIS 

We proceed to analyze the numerical performance of 
the semiclassical scattering factor (|66| as compared to 
the Glauber eikonal approximation ([57)1 and to the ele- 
mentary eikonal approximation used so far in evaluations 
of the Talbot-Lau effect. We will see that the elementary 
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25 30 35 40 

distance to the grating wall (nm) 



FIG. 5: Phase difFerence of the semiclassical (solid line) and 
the eikonal approximation (dotted line) to the exact result, 
plotted versus the distance to one grating slit wall. We use 
grating and particle parameters as in Fig. [2] (corresponding 
to the fuUerene experiment HC 



eikonal approximation was appropriate in the molecu- 
lar matter wave experiments performed to date [20, 122] ■ 
At the same time, both the Glauber and the semiclas- 
sical approximation significantly improve the treatment 
of laser gratings if the particles have larger polarizabil- 
ities or smaller velocities (as required if their mass ins 
increased). 

According to the general theory from Sect. IIIBI the 
interference pattern (P^ is determined by the general- 
ized, momentum dependent grating coefficients (j39p via 
the scattering factor ([55)1 . In order to assess the validity 
of the different approximations it is therefore pertinent 
to evaluate the scattering factors directly and compare 
them to a numerical implementation of the exact prop- 
agation of a plane wave through the grating interaction 
region. 



A. The scattering factor 



The exact transverse scattering factor (|38)) for the grat- 
ing can be obtained by a numerical evaluation of the limit 



(f) {x,p) = lim 



(a;|uL%U2TUL%|p) 



{x\p) 



(69) 



In practice, this is done by computing the propagator 
matrix elements by means of a split operator technique 
[49l |50 | , making sure that the propagation time T is suffi- 
ciently large (much larger than the grating passage time) 
so that the result is converged. In the following examples 
we use periodic boundary conditions for the position co- 
ordinate X and a fixed transverse momentum p not larger 
than the typical beam spread 10~^Pz. 



shows the phase differences between the approximate 
scattering factors and the exact scattering factor as a 
function of the distance to the wall. The calculations 
were done for a de Broglie wavelength of A = 4 pm. More- 
over, we choose an orthogonally incident plane wave, 
p = 0, so that the elementary and the Glauber eikonal 
approximations coincide. They are given by the dotted 
line, while the solid curve corresponds to the semiclassi- 
cal approximation (j66p . 

As one observes in Figure O the divergence of the wall 
potentials invalidates the approximations for the scatter- 
ing factor close to the walls. However, the contributions 
from this small vicinity of the slit walls do not affect the 
interference visibility appreciably, since it corresponds 
only to a small fraction of the semiclassical trajectories. 
This suggests that a reasonable cut-off criterion is given 
by the critical distance Xc to the wall, where a classical 
beam particle would hit the wall within the grating pas- 
sage time t — b/vz- Disregarding the initial transverse 
momentum p and the potential of the opposite wall, the 
time for a particle starting at the distance xq to hit the 
wall is given by 



T 



dx ■ 



1 



V^Vjx) 



(70) 



For a wall potential V{x) 
critical distance 



-C4X * this leads to the 



/ISmCib^ 



Pz 



1/6 



(71) 



The parameters used for Fig. [5] yield Xc — 21 nm, and 
this value indeed corresponds to the position where the 
semiclassical and the exact phase deviate by about 2tt. 
The fact that the eikonal approximation is virtually iden- 
tical to the exact phase factor for most of the slit width 
explains why the eikonal approximation is well justified 
with thin material gratings. In fact, our numerical re- 
sults indicate that in this case the eikonal approximation 
remains valid even for particles with a stronger particle- 
wall interaction, breaking down only in a regime where 
the interference visibility is already strongly diminished 
by the interaction effect. We therefore focus on the 
KDTLI setup in the following. 

As a final point, we note that the opening width of the 
slits is effectively reduced by twice the critical distance 
Xc- For large and slow particles it may be necessary to 
take this into account even at the first and at the third 
grating, since the fringe visibility depends quite sensi- 
tively on the corresponding effective open fractions. 



2. Laser grating 



1. Material grating 

We first take the grating to be of material type with a 
retarded Casimir-Polder interaction potential. Figure [5] 



Replacing the diffraction grating by a standing laser 
wave leads to the smooth and bounded interaction poten- 
tial ([M)l . For the numerical evaluation of the scattering 
factor, shown in Fig. [51 we choose the same parameters as 
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FIG. 6: Phase difference with respect to the exact scattering factor over one period d = 266 nm of the laser grating (parameters 
as in Fig.|3] with A = 3 pm, corresponding to the experiment [221). Panel (a) shows the results for an incoming plane wave with 
vanishing transverse momentum, p = 0, while (b) corresponds to p = 10~^ Pz. The solid line gives the error of the semiclassical 
approximation, which is not resolved on the scale given by the errors of the Glauber approximation (dashed line) and the 
elementary eikonal approximation (dotted line). The latter are indistinguishable in (a); notice the different scales of the y-axes. 



for Fig. [3] (motivated by the experiment [22]) and choose 
A = 3 pm. The longitudinal laser waist Wz = 20 /im is 
much larger than the grating period d = 266 nm. It fol- 
lows that non-zero transverse momenta p j^ now have 
to be considered separately since the free transverse mo- 
tion over the distance Wzp/pz must not be neglected. 

Figure [6] shows how the phases of the approximate 
scattering factors deviate from the exact phase. Panel 
(a) corresponds to a perpendicular incidence of the in- 
coming plane wave, p = 0, while (b) is evaluated for 
p = 10~^pz, as found in a beam spread of Imrad The 
Glauber and the elementary eikonal approximation coin- 
cide in the case (a) of perpendicular incidence, and they 
deviate from the exact phase by about 1 mrad. The er- 
ror of the semiclassical phase is smaller by three orders 
of magnitude, and it is not resolved in the plot. On the 
other hand, in the case (b) of a non-zero transverse mo- 
mentum the elementary eikonal approximation deviates 
substantially from the exact result, while the error for the 
Glauber approximation remains on the order of 10^"^ and 
the semiclassical one on the order of 10^^ (not resolved 
in the plot). 

As for the corresponding amplitude of the incident 
plane waves, the exact calculation yields deviations from 
the incident amplitude 1 on the order of 10^® in case (a) 
and deviations on the order of 10^^ in case (b), respec- 
tively (not shown). While the eikonal approximations 
cannot account for this effect, the semiclassical ampli- 
tude (|67p reproduces the exact result with an error of 
less than 10~®, and 10~^, respectively. 

The semiclassical expression of the scattering factor is 
thus demonstrated to be superior by orders of magni- 
tude compared to the eikonal approximations. While the 
Glauber eikonal approximation already improves the ele- 
mentary eikonal approximation significantly, it still does 
not take the amplitude modification into account. How- 
ever, the overall corrections to the eikonal approxima- 
tions are so small in the present parameter regime that 
the Talbot-Lau interference contrast is hardly affected. 



as demonstrated below. The elementary eikonal approx- 
imation thus remains valid for the considered experiment 

The situation changes distinctively if the de Broglie 
wavelength is increased by a factor of ten, to A = 30 pm. 
The resulting phase difference and amplitude plots are 
shown in Fig. [7] for the case of a transverse momentum 
p = 10~^pz. Now both the Glauber and the elementary 
eikonal phase strongly differ from the exact result, as 
demonstrated by the dashed and the dotted curves in 
Fig.EKa). At the same time, the semiclassical result (solid 
line) deviates by less than 100 mrad from the exact phase, 
and also the corresponding amplitude, seen in Fig.[7l^b), 
faithfully approximates the exact one. 

The semiclassical expression (p5|) starts to fail only if 
we decrease the beam velocity to such an extent that the 
trajectories get strongly deflected during the increased 
passage time. This is expected since the derivation 
assumes the corrections due to deflection to be small. 
Equation (j66p thus extends the eikonal approximation 
to the smaller beam velocities required for more massive 
particles, but it does not cover the whole semiclassical 
wavelength regime. 



B. The Talbot-Lau visibility 

We can now discuss how the improved treatment of 
the grating interaction effect affects the Talbot-Lau in- 
terference visibility. We focus again on the laser grat- 
ing setup demonstrated in [2^] . In the Figures [BHTUl the 
eikonal results are obtained by calculating the visibility 
P^ by means of the Talbot-Lau coefficients ([55)) . The 
semiclassical calculation implements the generalized for- 
mula for the interference pattern (P^ . where the semi- 
classical scattering factor (|66p enters by means of the 
generalized grating coefficients (15^1) . Note that, unlike 
in the elementary eikonal approximation, it is now es- 
sential to incorporate the angular distribution D (p/pz) 
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FIG. 7: (a) Phase diflerences with respect to the exact calculation like in Fig. |6l but for an increased de Broglie wavelength 
A = 30pm and p = 10~^Pz. The elementary eikonal approximation (dotted line) and the Glauber eikonal approximation 
(dashed line) now differ markedly from the exact phase result, while the semiclassical approximation (solid line) deviates by 
less that 100 mrad. (b) Corresponding amplitude modification as obtained from the semiclassical approximation (solid line) 
and the exact calculation (dashed line). 
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FIG. 8: Interference visibility as a function of the longitudinal 
waist Wz of the laser grating (for A = 3 pm, starting from the 
experimental value w^ — 20 /im 22]). The dotted line gives 
the result of the elementary eikonal approximation (|33|l , while 
the solid curve represents the semiclassical result (l66|l . 



FIG. 9: Fringe visibility as a function of the laser power Pl, 
starting from the experimental value Pl = 6W [22|. For 
considerably larger powers, or, equivalently, for larger polar- 
izabilities of the particles, the elementary eikonal approxima- 
tion (dotted line) deviates markedly from the semiclassical 
resuh (l66l). 



into the calculation. It is set here to be a Gaussian 
D{p/pz) oc exp (^—{p/pz)^ /2a^) with a realistic width 
of Of = 1 mrad. 

One immediate consequence of the dependence of the 
semiclassical and the Glauber approximations on the 
transverse momentum is demonstrated in Fig. [SJ where 
the interference visibility is plotted versus the longitudi- 
nal laser waist Wz, starting from the experimental value 
Wz = 20 /xm. The elementary eikonal approximation 
(dotted line) is independent of Wz due to the longitudinal 
integration in ([55)1 . The semiclassical result (solid line), 
which takes into account the transverse motion through 
the laser field, decreases with growing Wz- While the dif- 
ference between the approximations is negligible at the 
experimental value, it becomes significant for a larger 
laser focus. The Glauber approximation reproduces the 
semiclassical result up to a precision of 10~^ and is there- 
fore indistinguishable from the solid line Fig. [51 This 
implies that the visibility loss with growing waist is due 
to the free transverse motion through the laser grating. 



rather than due to a considerable deflection of the tra- 
jectories. 

A similar result is presented in Fig. [51 where we in- 
crease the laser power starting from the experimental 
value Pl = 6W. This is equivalent to increasing the 
polarizability of the particles, see ([M)) . One observes 
that the semiclassical result (solid curve) decreases more 
rapidly than the eikonal visibility (dotted curve) as the 
strength of the phase grating is increased. The Glauber 
approximation again matches the semiclassical curve. 

Finally, in Fig. [TUl we present the interference visibil- 
ity as a function of the de Broglie wavelength A around 
30 pm, corresponding to a tenfold smaller beam veloc- 
ity than in the experiment. In this case one observes 
at some wavelengths a considerable difference between 
the Glauber eikonal approximation (dotted line) and the 
semiclassical approximation (solid line). 

Note that the sharp dips, where the visibility drops to 
zero in Fig.llOland Fig. [51 indicate a shift of the interfer- 
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FIG. 10: Talbot-Lau interference visibilities for the laser grat- 
ing setup [22 ] with C70 fuUerenes. We use increased de Broglie 
wavelengths corresponding to velocities around v^ = 16 m/s 
[sH - The semiclassical approximation (solid line) yields val- 
ues which are clearly distinguishable from the Glauber eikonal 
approximation (dotted line) . The elementary eikonal approx- 
imation differs more strongly. 



ence pattern by half its period. One has to take this into 
account when averaging the visibility over the velocity 
distribution of the particle beam. For example, if this 
distribution was so broad that it amounts to averaging 
over two subsequent peaks in Fig. [TU] it would lead to 
almost zero visibility. 

We have seen that in general the validity of the 
eikonal approximation depends on the grating interaction 
strength, the passage time, the longitudinal velocity, and 
the distribution of the transverse momenta in the parti- 
cle beam. For the specific experiment [2^ the elementary 
eikonal approximation breaks down by either decreasing 
the longitudinal velocity by a factor of 10, or similarly 
by increasing the laser power Pl, the longitudinal waist 
Wz, or the particle polarizability by an order of magni- 
tude. The latter are effects mainly due to the transverse 
motion of the beam particles, which is not taken into 
account by the elementary eikonal approximation. Here 
the Glauber eikonal approximation (|57p . i.e., the semi- 
classical scattering factor without deflection, would be 
already sufficient, at least in the experimentally accessi- 
ble regime. This is no longer the case if one increases the 
wavelength, since the deflection effect is more sensitive 
to the wavelength than to the grating parameters, limit- 
ing the validity of both the elementary and the Glauber 
eikonal approximation. 



VI. CONCLUSIONS 

We presented a general theory of the coherent Talbot- 
Lau interference effect. It allows to incorporate the in- 
teraction between particle and grating structure-a dom- 
inant effect in near field interference-at various degrees 
of approximation. Our treatment shows that it is neces- 
sary to account for detailed beam characteristics, such as 



the angular distribution, whenever one is required to go 
beyond the elementary eikonal approximation, or if one 
wants to quantify the experimental adjustment require- 
ments. 

Using the phase space formulation of quantum me- 
chanics, we identify the appropriate generalization of the 
Talbot-Lau coefficients. They serve to incorporate the 
most general coherent grating transformation and to de- 
scribe the various near field interference effects in a trans- 
parent fashion. The general effect of the passage through 
a grating can thus be formulated in terms of scattering 
theory, providing a starting point for the numerically ex- 
act evaluation of the interference pattern. 

Moreover, the semiclassical approximation of the S- 
matrix yields a systematic and non-perturbative im- 
provement over the elementary eikonal approximation. 
An additional high-energy approximation of the semiclas- 
sical trajectories then yields the Glauber eikonal approx- 
imation and the semiclassical deflection approximation 
as systematic corrections to the standard treatment. A 
comparison with the numerically exact calculation ver- 
ifies the high quality of the semiclassical deflection ap- 
proximation. It suggests that a Kapitza-Dirac Talbot- 
Lau interferometer, where the center grating is replaced 
by a standing light wave, will be able to demonstrate the 
wave nature even of particles which are so large that the 
eikonal approximation is no longer valid. 
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APPENDIX: CLASSICAL DESCRIPTION 

If the Talbot-Lau experiment is to prove the quantum 
nature of particles one clearly needs to be able to distin- 
guish between the quantum interference effect and the 
moire-type shadow effect that may occur with classical 
particles. A formulation is therefore required that yields 
the classical shadow contrast by using the same assump- 
tions and approximations as in the quantum case. We 
present this classical theory in the following by assuming 
a material grating of thickness b with a transverse inter- 
action potential V{x) of the grating slit walls. The case of 
an explicitly ^-dependent interaction potential V {x,z), 
such as a laser grating [22|, can be treated in the one- 
dimensional model by an effectively time-dependent po- 
tential V {x,t) = V {x,Pzt/m), where the longitudinal 
motion provides the time coordinate for a given momen- 
tum Pz ■ 

The classical formulation is based on the phase space 
density rather than the Wigner function. The corre- 
sponding classical propagator through a diffraction grat- 
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ing is given by the expression 



Kci{x,p;xq,po) = \t{xo)\ S [xq- Xg [x,p. 
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x,p, 



(A.l) 



where the hard grating wall cutoff of the particle beam is 
taken into account at the entrance into the grating. The 
phase space coordinates xf,p'^ are the starting point of 



the classical trajectory evolving to x,p under the influ- 
ence of the interaction potential V{x) within the grating 
passage time t = mb/pz. 

Since the free evolution of the phase space density is 
given by the same transformation ^ as in the case of the 
Wigner function, one ends up with the general classical 
shadow pattern, denoted by /a instead of w^, , after sub- 
stituting the classical propagator (jA.ip into the general 
Talbot-Lau calculation from Sect. Ill 51 
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One can numerically implement this formula directly, 
rather than performing a Fourier decomposition with re- 
spect to the argument x—prfL/pz of the trajectory terms. 
The ideal moire shadow pattern is obtained if one dis- 
regards both the interaction potential and the grating 
thickness by setting Xq' {x,p, t) — x and p'^ {x,p, t) — p. 
If a particle- wall interaction is present, the classical ana- 
logue to the eikonal approximation (j33|) is to approximate 
the deflection of the trajectory due to the grating by an 
instantaneous momentum kick 12611 . 
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The classical momentum kick coefficients read as 
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(A-5) 
If the second grating is implemented by a standing 
laser beam the classical calculation yields an analyti- 
cal expression for the Talbot-Lau coefficients B"^ (^) = 
J2k^m-kiO) Ck{£,)- They are related to the quantum 



expression (I35p by replacing the sine function in the ar- 
gument with its linear expansion, 



Putting this into the classical formula (jA.2[) . performing 
all the Fourier decompositions, and focusing on a partic- 
ular r : s Talbot-Lau resonance, as done in the eikonal 
quantum case ((28)) , one obtains the classical shadow pat- 
tern in eikonal approximation. 
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Since the argument of the _B,„ is proportional to the de 
Broglie wavelength A in the Talbot-Lau interference ef- 
fect, this means that the quantum interference and the 
classical shadow effect become indistinguishable in the 
naive classical limit of a vanishing wavelength, A — > 0. 
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